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ABSTRACT 

An  analytical/computational  technique  previously  developed  for  determining 
the  geometrically  and  constitutively  nonlinear  response  of  a  submerged, 
infinite  cylindrical  shell  to  a  transverse,  transient  acoustic  wave  is  used 
to  study  the  damage  behavior  of  the  shell.  Incident  waves  of  rectangular 
pressure-profile  are  considered,  nonlinear  transient  response  computations 
are  performed,  and  damage  results  are  described  in  terms  of  iso-damage 
curves  based  on  extensional  set  strain.  Results  generated  through  the  use 
of  the  doubly  asymptotic  approximation  for  treatment  of  the  fluid-structure 
interaction  differ  appreciably  from  their  exact  counterparts. 
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Section  1 
INTRODUCTION 

An  analytical/computational  technique  has  recently  been  developed  for  deter¬ 
mining  the  geometrically  and  constitutively  nonlinear  response  of  an  infinite 
cylindrical  shell  to  a  transverse,  transient  acoustic  wave  [1].  The  technique 
involves  the  use  of  the  structural  analyzer  DYNAPLAS  [2]  to  treat  shell  re¬ 
sponse,  and  the  use  of  the  residual  potential  formulation  [3]  to  treat  the 
fluid-structure  interaction;  it  has  been  implemented  in  the  form  of  a  coupled 
software  assembly  named  RPF-DYNA. 

In  this  study,  RPF-DYNA  has  been  used  to  investigate  the  damage  behavior  of 
a  particular  shell  for  incident  plane  waves  of  rectangular  pressure-profile. 

The  shell  exhibits  elastic/perfectly  plastic  material  behavior  and  is  charac¬ 
terized  by  hydrostatic  elastic-critical-buckling  and  elastic-limit  pressures 
that  are  virtually  equal.  Because  the  deformational  response  of  the  shell 
is  dominated  by  extensional  motion,  damage  results  are  described  in  terms  of 
iso-damage  curves  [4]  based  on  extensional  set  strain. 

Iso-damage  curves  are  also  constructed  for  the  shell  that  pertain  to  the  use 
of  the  doubly  asymptotic  approximation  (DAA)  [5,6]  for  treatment  of  the  fluid- 
structure  interaction.  This  approximation,  which  is  the  basis  for  fluid- 
structure  interaction  analysis  in  a  number  of  existing  codes  [7-12],  is 
asymptotically  exact  for  both  low-  and  high-frequency  fluid  motions,  effect¬ 
ing  a  smooth  transition  in  the  intermediate  frequency  range.  Its  computational 
advantage  is  that  it  may  be  expressed  as  a  matrix  ordinary  differential  equa¬ 
tion  without  requiring  discretization  of  the  infinite  volume  of  fluid  sur¬ 
rounding  the  structure. 
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Section  2 

RESPONSE  CALCULATIONS 


Consider  the  two-dimensional,  plane-strain  motions  of  the  submerged,  infinite, 
circular  cylindrical  shell  of  Figure  1.  The  shell  is  excited  by  a  transient 
acoustic  wave  that  first  contacts  the  shell  at  e  =  n..  During  the  resulting 
fluid-structure  interaction,  shell  behavior  may  involve  both  geometric  and 
constitutive  nonl inearity. 


2.1  NONDIMENSIONAL  RESPONSE  EQUATIONS 


Following  an  expansion  of  the  pertinent  response  variables  in  circumferential 
Fourier  series,  nondimensional  response  equations  for  each  Fourier  harmonic 
may  be  derived  as  [1] 


(VH^o  +  *  +  fo(v’w)  =  *Io  +  uIo  '  \  *$o 
(y)(?)Wn  +*n  +  fwn(v’w)  =  *In  +  uIn  “  I  * 

+  (ilo  +  *So)(nv 


+  *Ro 
Sn  +  *Rn 


n  >  1 


(1) 


(f)(i)Vn  +  fvn(v’w)  =  (*Io  +  *So)(vn  +  nwn}’  n  ^  1 

*Sn  +  1  *Sn  =  uIn  '  wn  +  *Rn 
*Rn  "  "rn  ^Sn 


where  v  and  w  are  circumferential  and  radial  shell -displacement  harmonics, 
n  n 

respectively,  f0»  fvn  and  fwn  are  stiffness-force  harmonics  computed  within 
DYNAPLAS  that  involve  1  inear-elastic,  geometrically  nonlinear  and  constitutively  non¬ 
linear  behavior,  and  UjN  are  fluid-velocity-potential  and  radial-fluid-particle- 
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velocity  harmonics  for  the  incident-wave  field  at  the  shell's  wet  surface, 
respectively,  and  <pRn  are  wet-surface  scattered-wave-velocity-potential  and 
residual-velocity-potential  harmonics , respectively ,  rp  is  a  known  time- 
dependent  characteristic  function  [3],  and  where  n  is  the  Fourier  index,  a 
dot  denotes  temporal  differentiation  and  the  asterisk  denotes  temporal  con¬ 
volution.  All  of  the  quantities  and  operations  in  (1)  are  nondimensional , 
normalized  in  accordance  with  the  convention 

Vd.  •  w/a'  Vd.  “  ct/a-  Vd.  =  */ac>  Vd.  *  f/cc2  <2> 

For  n=0,  (1)  contains  three  equations  for  the  three  unknowns  wQ,  4>Sq  and 
^Ro’  w^e  ^or  e‘>c^  n  >  1*  (1)  contains  four  equations  for  the  four  unknowns 
vn’  wn’  ^Sn  and  ^Rn'  Here’  as  in  the  first  three  of  (1)  are  solved 
with  the  half-step  central -difference  algorithm,  the  fourth  of  (1)  is  solved 
with  a  fourth-order  Runge-Kutta  scheme,  and  the  last  of  (1)  is  solved  by 
trapezoidal  integration.  This  procedure  possesses  satisfactory  accuracy, 
stability  and  efficiency  characteristics. 

2.2  DAMAGE  BEHAVIOR 

The  damage  behavior  of  a  particular  cylindrical  shell  with  elastic/perfectly 
plastic  material  behavior  is  studied  by  calculating  numerous  transient  re¬ 
sponse  histories  for  incident  plane  waves  of  rectangular  pressure-profile. 

The  shell  chosen  is  the  single-layer  "compromise  shell"  of  [1],  whose  inertial, 
elastic,  hydrostatic-elastic-stability,  and  extensional  elastic-limit  charac¬ 
teristics  match  those  of  a  steel  sandwich  shell,  which  exhibits  in  plane- 
strain  fashion  the  enhanced  flexural  stiffness  properties  of  stiffened  shells. 

As  discussed  in  [1],  the  use  of  the  compromise  shell  is  necessitated  by  the 
limitation  of  DYNAPLAS  analysis  capability  to  single-layer  shells.  Table  1 
shows  the  basis  of  equivalence  for  the  two  shells;  the  only  discrepancies 
pertain  to  inelastic  flexural  characteristics,  which,  as  seen  below,  are  of 
minor  significance.  Note  that  the  yeld-stress  value  ay  is  such  that  the  shells' 
hydrostatic  elastic-limit  pressure  P0  (whose  non-dimensional  value  is  identical 
to  the  elastic-limit  membrane  stress-resultant  Ny)  exceeds  their  static-elastic- 
stability  pressure  P0  by  2.4%. 
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The  response  variable  chosen  for  damage  assessment  is  circumferential  exten- 
sional  (or  membrane)  set  strain,  for  two  reasons.  First,  as  discussed  in  [1], 
the  inelastic  response  of  this  shell  is  dominated  by  extensional  motion.  Second, 
extensional  shell  response  almost  always  reaches  its  late-time  asymptotic 
limit  within  five  shell -envelopment  periods  following  the  passage  of  the  in¬ 
cident  wave  over  the  shell;  this  is  in  contrast  to  flexural  shell  response, 
which  exhibits  low-frequency  oscillatory  behavior  for  extremely  long  times. 

Hence  a  systematic  treatment  of  flexural  set  strain  requires  that  each  tran¬ 
sient  response  calculation  be  followed  by  a  dynamic-relaxation  calculation  in¬ 
volving  the  introduction  of  artificial  damping  to  damp  out  the  oscillations. 

In  view  of  the  first  reason,  the  results  thus  obtained  are  probably  not  worth 
the  effort. 

Representative  inner-fiber  and  outer-fiber  strain  response  histories  are 

shown  in  Figure  2.  In  this  case,  the  magnitude  of  the  incident  rectangular 

wave  is  five  times  the  static-elastic-stability  pressure  (Pj  =  5  P  )  and  the 

duration  of  the  wave  is  equal  to  the  shell -envelopment  period  (Tj  =  2).  For 

these  input  parameters,  extensional  set  strain  is  largest  in  magnitude  at 

9=0.  From  Figure  2,  |e^e^|  =  4.8%,  while  |e«*  |  ~  0.15%;  hence  peak 

2Xl  max  t  i  ex 

flexural  strain  is  only  3 %  of  peak  extensional  strain  in  this  case. 

Highly  unrepresentative  inner-fiber  and  outer-fiber  strain  response  histories 
are  shown  in  Figure  3.  In  this  case,  for  which  P j  =  1.8  Pc  and  Tj  =  10,  ex¬ 
tensional  set  strain  is  largest  in  magnitude  at  9  =  708.  From  the  figure, 

|e^e^|  =  6.0%,  while  I  «  2.5%;  hence  peak  flexural  strain  is  about 

ext  max  T  lex  rpax 

40%  of  peak  extensional  strain  in  this  case. 

A  more  comprehensive  picture  of  the  inelastic  response  results  obtained  in 
this  study  is  provided  in  Table  2.  This  table  summarizes  the  results  obtained 
in  18  of  the  more  than  50  RPF-DYNA  transient  response  calculations  performed. 

Six  quantities  are  shown  for  each  response  calculation,  as  indicated  in  the 
format  statement. 

Table  2  shows  that  the  set-strain  field  in  the  shell  is  dominated  by  the  n  =  0 
harmonic,  with  the  n  =  1  and  n  =  2  extensional  harmonics  making  modest 
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contributions  and  the  higher  extensional  harmonics  contributing  very  little. 
Flexural  set  strain  is  significant  only  in  a  few  cases.  Maximum  extensional 
set  strain  usually  occurs  on  the  side  of  the  shell  that  faces  away  from  the 
incoming  wave. 

2.3  I SO- DAMAGE  CURVES 

Calculated  values  of  le^|  are  shown  as  circled  dots  in  Figure  4,  plotted 

ext  ma,< 

as  functions  of  magnitude  (Pj  and  duration  (Tj)  of  the  incident  rectangular 

wave.  (Recall  that  Pc  isthe  elastic-critical-buckl ing  pressure  for  this  shell, 

which  is  nearly  equal  to  its  hydrostatic  elastic-limit  pressure).  The  Tj  = 

const,  curves  that  connect  the  dots  permit  the  accurate  estimation  of  six 

combinations  of  PT  and  Tt  that  yield  a  prescribed  value  of  |e  +1  .  The 

II  ext 'max 

selection  of  five  such  values  then  leads  to  the  iso-damage  curves  of  Figure  5. 

Figure  5  shows  curves  of  constant  plotted  in  terms  of  the  pressure 

magnitude  and  total  impulse  that  completely  characterize  an  incident  rectangu¬ 
lar  wave.  Although  the  curves  exhibit  the  general  characteristics  of  iso¬ 
damage  curves,  i.e.,  transition  from  vertical  asymptotes  for  short  pulses  to 
horizontal  asymptotes  for  long  pulses,  they  do  not  possess  the  simple  hyber- 
bolic  shapes  produced  by  simple  mechanical  systems  [4], 

In  Figure  5,  each  iso-damage  curve  appears  to  consist  of  two  damage  branches, 
one  for  pulses  with  Tj  <"  3  and  another  for  pulses  with  Tj  >  3.  The  distinctly 
positive  slopes  characterizing  the  high-strain  Tj  <  3  branches  for  Tj  >  1.5 
produce  an  unexpected  characteristic,  viz.,  that  a  decrease  in  pressure  magni¬ 
tude,  PT,  for  fixed  impulse,  P.  Tt,  may  produce  an  increase  in  peak  extensional 

^  .  i  S6t |  *  ^ 

set  strain,  |e  J  .  This  is  undoubtedly  associated  with  a  resonance  condi- 

ext  max 

tion  where  the  shell  is  especially  susceptible  to  pulses  with  durations  com¬ 
parable  to  the  shell -envelopment  period. 

2.4  DOUBLY  ASYMPTOTIC  APPROXIMATION 

As  discussed  earlier,  the  doubly  asymptotic  approximation  (DAA)  has  been  exten¬ 
sively  used  in  problems  of  this  type  for  treatment  of  the  fluid-structure 
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Extensive  DAA-DYNA  calculations  have  been  performed  for  the  compromise  shell 
of  Table  1  The  results  of  these  calculations  exhibit  dominance  of  the  n  =  0 
harmonic,  as  previously  seen  in  the  RPF-DYNA  results;  however,  n  =  1  strain  re¬ 
sponse  emerges  in  the  DAA-DYNA  calculations  as  almost  as  significant  as  n  =  0 
strain  response.  Contributions  to  total  strain  by  the  higher  harmonics,  whether 
extensional  or  flexural,  are  generally  small  in  the  DAA-DYNA  calculations. 

Values  of  |efe^|  obtained  from  the  DAA-DYNA  calculations  are  shown  as  circled 
exi-  max 

dots  in  Figure  6,  plotted  in  the  manner  of  Figure  4.  In  contrast  to  their 
RPF-DYNA  counterparts,  these  results  indicate  that  peak  extensional  set  strain 
usually  occurs  on  the  side  of  the  shell  that  faces  toward  the  incoming  wave. 


The  same  procedure  that  produced  the  iso-damage  curves  of  Figure  5  from  the 
Tj  =  const,  curves  of  Figure  4  has  been  used  to  produce  the  iso-damage  curves  of 
Figure  7  from  the  Tj  = const,  curves  of  Figure  6.  Even  a  cursory  comparison  of 
Figures  5  and  7  reveals  marked  differences  between  the  RPF-DYNA  and  DAA-DYNA 
curves.  First,  and  most  important,  the  DAA-DYNA  curves  for  higher-strain  values 
lie  well  to  the  right  of  and  above  their  RPF-DYNA  counterparts,  which  means  that 
use  of  the  DAA  may  lead  to  serious  underpredictions  of  damage  levels.  Second, 
and  less  important,  the  DAA-DYNA  curves  do  not  exhibit  the  resonance  behavior 
observed  in  the  RPF-DYNA  curves. 


Figs.  5  and  7  demonstrate  the  general  inadequacy  of  the  DAA  for  treatment  of 
the  fluid-structure  interaction  in  the  present  problem.  As  mentioned  in  [1], 
the  reason  for  this  inadequacy  is  that  the  DAA  reduces  to  the  PWA  for  n  =  0 
response,  which  severely  attenuates  the  rather  low-frequency  axisymmetric 
motion  characterizing  the  inelastic  response. 

To  gain  insight  into  possible  improvement  of  the  situation  just  described,  re¬ 
sponse  calculations  for  Pj  =  4Pc,Tj  =  4  were  performed  using  the  cylindrical 
wave  approximation  (CWA)  [13].  Utilization  of  this  approximation  in  the  present 
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problem  merely  involves  taking  ^  ^*0in(l).  Displacement  response  histories 
for  n  =  0  and  n  =  1  extensional  motion  produced  by  RPF-DYNA,  DAA-DYNA  and 
CWA-DYNA  are  shown  in  Figure  8.  It  is  seen  that  the  RPF-DYNA  histories  are 
bounded  by  their  DAA-DYNA  and  CWA-DYNA  counterparts;  unfortunately,  the  bounds 
are  unacceptably  large.  Similar  histories  for  n  =  1  and  n  =  2  flexural  motion 
are  shown  in  Figure  9.  Here,  the  RPF-DYNA  histories  are  bounded  by  their 
DAA-DYNA  and  CWA-DYNA  counterparts  over  most  of  the  time  span  shown;  these 


bounds  too  are  unacceptably  large. 
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Section  3 
CONCLUSION 


The  extensive  calculations  performed  in  this  study  have  provided  a  rather  com¬ 
prehensive  picture  of  the  damage  behavior  of  a  submerged,  infinite  cylindrical 
shell  excited  by  transverse  acoustic  waves  of  rectangular  pressure  profile.  This 
behavior  is  summarized  in  the  iso-damage  curves  of  Figure  5. 

The  calculations  have  also  confirmed  the  failure  of  the  doubly  asymptotic  ap¬ 
proximation  as  a  satisfactory  treatment  of  the  fluid-structure  interaction  for 
the  inelastically  excited,  infinite  cylindrical  shell.  Such  failure  was  first 
indicated  by  a  sigle  response  calculation  reported  in  [1].  A  similar  calculation 
based  on  the  cylindrical  wave  approximation  has  produced,  in  this  study,  slightly 
better,  but  still  unsatisfactory,  results. 

The  reason  for  the  failure  of  the  doubly-asymptotic  approximation  is  simply 
that,  for  the  axisymmetric  response  of  the  infinite  cylindrical  shell,  it  reduces 
to  the  plane  wave  approximation,  which  is  only  singly  asymptotic.  Such  reduction 
does  not  occur,  however,  for  finite  three-dimensional  bodies;  in  fact,  the  OAA 
is  exact  for  the  n  =  0  response  of  a  spherical  body. 

On  the  assumption  (which  has  yet  to  be  verified)  that  the  DAA  is  satisfactory 
for  determination  of  the  inelastic  response  of  a  submerged  spherical  shell  to 
transient  acoustic  waves,  the  following  question  arises:  How  large  may  the 
aspect  (length/diameter)  ratio  become  before  the  DAA  fails  as  a  satisfactory 
treatment  of  the  fluid-structure  interaction  for  inelastic  shell -response 
problems?  Seeking  the  answer  to  this  question  is  a  fitting  subject  of 
future  work. 
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Table  L.  SINGLE-LAYER  COMPROMISE  AND  SANDWICH  SHELLS 


Single-Layer  Shell 


Sandwich  Shell 


IT 

%  h 

11 


H  Middle  Layer  Contributes 

i  No  Mass  or  Stiffness 

_L_  77  /  Trrrr  nrrri  riTrrrrrrrri 


-  -  0.01,  -  -  0.06266 
a  a 


—  =  0.772,  —  *  3.53,  — -  0.02188  — -  7.72,  —  -  3.53,  — *  0.2188 

P  c  ’2  p  c  2 

pc  K  pc 


Inertial  and  Elastic  Properties 


2  2 
o  h  P  C  h  Dnh  P  C  h 

":t" 0.0772,  "i-t"  =  *  0.9620  "pt" 0.0772,  ”Et”  -  — ■ y-  *  0.9620 

pa  *  2  pa  2 

pc  a  pc  a 


"El"  *  •  y^(-)  -  8.017xl0"4 

l  u  a 

CC 


„-T„  -  ^0  .  100  ,h 3 


f  ’  W?  “  8-017xl° 


Static-Elastic-Stability  and  Exfensional  Elastic-Limit  Characteristics 


a  c  2  ,  3  . 

1  *  k(THc)  (I}  *  2-405*10 


p  c  2  .3  , 

Pr  -  25(~)(-f)  e)  -  2.4)5 xl0“J 
c  oca 


,  o  ~ 

(7U-JI)  *  ^.^ixio--1  ** 


,o  . 

N  -  (-)(— 4)  *  2 . 462x10“ J  ** 
y  a  2. 

J  pc 


Inelastic  Flexural  Characteristics 


My  -  4.  103x10 
Ky  -  0.05113 
Mu  -  6.155x10“ 


Yield  Moment 
Yield  Curvature 
Ultimate  Moment 


M  *  6.549x10 

y 

\  ■  0.08170 

M  -  7.099x10“ 
u 


See  Figure  1  for  parameter 
definitions 


von  Mises  yield  condition 
with  v  *  0.3 


Table  2.  RPF-DYNA  RUNS:  SET  STRAIN  SUMMARY 


Format 

Maximum  Extenslonal  Set  Strain  Magnitude,  |e**i|  P,  ■  P,/Pp 

*xl  max  1  *  c 

Circumferential  Location  Where  Maximum  Occurs 

Percentage  Contribution  of  n«0  Set  Strain  to  |e**J| 

ext  max  ,  set, 

Percentage  Contribution  of  n*l  Extenslonal  Set  Strain  to  |e_„J 

eXt 

Percentage  Contribution  of  n«2  Extenslonal  Set  Strain  to  |e|f|| 

,  set,  ■  set,  exl  m*x 

Approximate  Ratio  of  to 


TI  *  *•  PI  * 

12 

V 

i.  \  *  20 

A 

PI 

*  28 

1.823% 

5.659% 

11. 39% 

120° 

0° 

0° 

78% 

66% 

61% 

11% 

9% 

19% 

12% 

16% 

15% 

0.06 

0.02 

0.00 

A 

A 

A 

V'-V 

3.6 

V 

1,  P,  *  12 

t_l 

■l.P, 

*  18 

1.089% 

6.159% 

12.40% 

0° 

0° 

70° 

41% 

78% 

87% 

16% 

10% 

5% 

20% 

9% 

7% 

0.15 

0.01 

0.27 

A 

A 

^  A 

TI  *  2’  PI  * 

3 

V 

2.  P,  *  5 

II 

PJ 

*  9 

1.861% 

4.801% 

12.09% 

0° 

0° 

0° 

55% 

63% 

80% 

24% 

14% 

18% 

12% 

13% 

0% 

0.06 

0.03 

0.06 

A 

A 

\ 

TI  “  4’  PI  * 

1.8 

V 

4,  Pj  *  3.6 

t_l 

*  4,  Pj 

•  4.2 

1.388% 

8.305% 

11.09% 

60° 

0° 

0° 

69% 

64% 

66% 

7% 

23% 

23% 

8% 

10% 

9% 

0.07 

0.02 

0.01 

A 

A 

A 

Ti  * 6-  V 

1.7 

V 

6,  Pj  «  2.2 

Ll 

•6.  P, 

•  3 

1.890% 

4.730% 

11.09% 

60° 

50° 

0° 

89% 

73% 

66% 

0% 

20% 

27% 

5% 

2% 

7% 

0.04 

0.08 

0.00 

A 

A 

A 

TI  »  10,  Pt  ■ 

■  1.2 

V 

10,  Pj  *  1.8 

II 

*  10,  P, 

,-2.2 

1.182% 

5.972% 

12.41% 

180° 

70° 

40° 

83% 

89% 

78% 

11% 

3% 

20% 

5% 

8% 

2% 

0.10 

0.42 

0.41 
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EXTREME-FIBER  STRAINS,  e(0,  +z 


0.0  1.0  8.0  12.0  16.0  20.0  21.0  28.0 


NONDIMENSIONAL  TIME,  t 


Figure  2  Representative  RPF-DYNA  Strain  Response  of  Compromise 
Shell  (IF  denotes  inner  fiber;  OF  denotes  outer  fiber; 
?z  «  4PC  4P0  ,  T1  -  4) 
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EXTREME-FIBER  STRAINS,  e(6, 


EPS [60 ) / I F 
EPS (60) /OF 
EPS (70) /IP 
EPS (70) /OF 
EPS (80) /IF 
EPS (80) /OF 


Figure 


NONDIMENS IONAL  TIME,  t 

3  Unrepresentative  RPF-DYNA  Strain  Response  of  Compromise 
Shell  (IF  denotes  inner  fiber;  OF  denotes  outer  fiber, 

P 


2  4  6  8  10  20 

VPc 


Maximum  Extens ional  Set- Strain  Magnitudes  as 
Calculated  by  RPF-DYNA  (angles  indicate  the 
locations  at  which  the  maxima  are  reached) 
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pITI/P 
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Figure  5  RPF  Iso-Damage  Curves  from  the  Data  of  Figure  4 
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Figure  7  DAA  Iso-Damege  Curves  from  Che  Dace  of  Figure  6 
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NORMALIZED  DISPLACEMENT,  -w„(t)/PT  AND  e,  (t)/P 


Figure  8  n=0  and  n«l  Extensional  Displacement  Response  of 

Compromise  Shell  as  Computed  with  RPF-DYNA,  DAA-DYNA 
and  CWA-DYNA  (Pj  -  4PC  »  4P0  ,  Tj  -  4) 


20 


NORMALIZED  DISPLACEMENT,  f, (t)/PT  AND  -f9(t)/P 


Figure  9  n~l  and  n*2  Flexural  Displacement  Response  of  Compro¬ 
mise  Shell  as  Computed  with  RPF-DYNA,  DAA-DYNA  and 
CWA-DYNA  (Pj  .  APC  *  4P0  ,  T].  =  4) 
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